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Abstract 

The magnetic structure and fluctuations of tetragonal GdRhlns were studied by resonant x-ray- 
diffraction at the Gd Ljj and Lm edges, followed by a renormalization group analysis for this and 
other related Gd-based compounds, namely Gd2lrlns and Gdln3. These compounds are spin-only 
analogs of the isostructural Ce-based heavy-fermion superconductors. The ground state of GdRhlns 
shows a commensurate antiferromagnetic spin structure with propagation vector f = (0, \,\), 
corresponding to a parallel spin alignment along the a-direction and antiparallel alignment along b 
and c. The spin direction lies along a. A comparison between this magnetic structure and those of 
other members of the i? m (Co,Rh,Ir) n In3 m +2n family (R =rare earth, n = 0, 1; m = 1,2) indicates 
that, in general, f is determined by a competition between first- (Ji) and second-neighbor (J2) 
antiferromagnetic (AFM) interactions. While a large J1/J2 ratio favors an antiparallel alignment 
along the three directions (the so-called G-AFM structure), a smaller ratio favors the magnetic 
structure of GdRhlns (C-AFM). In particular, it is inferred that the heavy-fermion superconductor 
CeRhlns is in a frontier between these two ground states, which may explain its non-collinear 
spiral magnetic structure. The critical behavior of GdRhlns close to the paramagnetic transition 
at T/v = 39 K was also studied in detail. A typical second-order transition with the ordered 
magnetization critical parameter (5 = 0.35 was experimentally found, and theoretically investigated 
by means of a renormalization group analysis. Although the Gd 4/ 7 electrons define a half-filled, 
spherically symmetrical shell, leading to a nearly isotropic spin system, it is argued that a significant 
spin anisotropy must be claimed to understand the second order of the paramagnetic transition of 
GdRhlns and the related compound Gd2lrlns. 

PACS numbers: 75.25,+z, 75.40. Cx, 75.50.Ee, 61.10.Nz 
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I. INTRODUCTION 



The recent discovery of a new class of heavy-fermion superconductors, 
Ce m (Co,Rh,Ir) n In 3m+2n (n = 0,1; m = 1, 2) , 1 i 2 i 3 i 4 i 5 has triggered extensive research 
on their physical properties. While it is well established that the spin fluctuations are key 
ingredients to determine the superconductivity in this family— the details of the pairing 
mechanism are not presently known.— Nonetheless, the relatively large critical temperatures 
found for some compounds (T c = 2.3 K for tetragonal CeCoIns is a record-high value for 
Ce-based heavy fermion systems), and the evolution of this property with the ratio of c 
and a lattice parameters^ indicates that ingredients closely related with the crystalline 
environment must be taken into account. This connection can be further explored by 
direct investigations of the sensitivity of the electronic structure on small deviations of the 
crystalline environment, as well as of the influence of the such environment on the nature 
and magnitude of the spin structures and the fluctuations that presumably mediate the 
superconductivity. In line with the second approach, the magnetic structures below T/v 
were resolved for a number of compounds, revealing interesting trends. The cubic Cel ^ 10 ' 11 
and tetragonal Ce2RhIngJ^ show commensurate magnetic structures with antiferromagnetic 
(AFM) alignment of Ce spins along the three nearest-neighbor directions (here called 
G-AFM phase in analogy to the nomenclature used in the manganites 13 -), and CeRhlns 
forms an incommensurate spiral along the c-direction while still keeping the AFM coupling 
in the a&-plane . Ml15 Recently, coexisting magnetic orders were observed in the alloy system 
CeKh.i_ x lY x ln 5 in the interval 0.25 < x < 0.6 where AFM and superconducting phases 
overlap.— The first AFM phase is identical to the spiral phase of CeRhIn 5 , while the second 
one shows AFM alignment along the three nearest-neighbor directions such as in Celn a 10 ' 11 
and Ce 2 Riling.— These results indicate that distinct AFM ground states compete in this 
system, and that such competition may favor or be related with superconductivity. 

Magnetic phenomena and their connection to crystal structures can be further explored 
by a thorough investigation of other members of the i? m (Co,Rh,Ir) n In3 m+ 2 n family (R =rare 
earth ^ Ce, n = 0,1; m = 1,2). Since such compounds are not heavy fermions and/or 
superconductors, the knowledge thus obtained may be taken as a starting point to under- 
stand the more complex and rich behavior of the Ce-based compounds. Previous macro- 
scopic studies on the i? m (Co,Rh,Ir)„In 3m+2n family (R = Pr, Nd, Sm, Gd, and Tb) indicate 
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that the evolution of the Neel temperature (Tjv) with R does not follow the de Gennes 
scaling , 17 ! 18 suggesting that crystal field or other anisotropy effects may be important to 
determine the critical temperatures and perhaps even the magnetic ground state in this sys- 
tem. Neutron scattering studies on Ndln 3 indicate a ground state with AFM alignment along 
two nearest-neighbor directions and FM alignment along the third direction, the so-called 
C-AFM structure.— Additional phases with modulated moments along the FM directions 
have also been identified below the Neel temperature for this compound . 19 i 20 Similar C-AFM 
ground states have been identified in NdRhlns 21 , Gd2irln s 22 and TbRhIn 5 .— 

In continuation of the attempt to build a minimum comprehension of the magnetism of 
a spin-only system (L = 0) under similar crystal environments to Ce m (Co,Rh,Ir) n In3 m+2n 
(n = 0, 1; m = 1, 2) heavy- fermion superconductors, we carried out an experimental and the- 
oretical investigation of the magnetic structure and fluctuations of a GdRhlns single crystal. 
We note that the magnetism of Gd-based compounds is expected to be particularly simple, 
due to the absence of large orbital moments and crystal-field interactions associated with a 
half filled 4/ 7 shell. The magnetic structure resolved here is of the C-AFM type, with partly 
frustrated first-neighbor (Ji) spin interactions. Based on this result and previous studies on 
other members of the i? m (Co,Rh,Ir) n In 3m+ 2 n family, we conclude that the relative orienta- 
tion between neighboring i?-spins is determined primarily from a close competition between 
first- (Ji) and second-nearest-neighbor (J2) AFM exchange interactions. This may be an 
important factor behing the complex magnetic behavior of the Ce-based compounds. The 
critical behavior close to the magnetic ordering transition in GdRhlns was also investigated. 
The magnetic order parameter shows a power-law behavior close to the Neel temperature, 
characteristic of a second-order transition such as in Gd2lrln 8 . 22 A theoretical renormaliza- 
tion group analysis was performed. Our results suggest that spin anisotropy terms, possibly 
arising from dipolar and other fairly weak interactions, must be claimed to understand the 
second order of the paramagnetic transition of GdRhIn 5 and the related compound Gd2lrln 8 . 

II. EXPERIMENTAL DETAILS 

A single crystal of GdRhlns was grown by the In-flux method as described previousl y 17124 
The studied surface was finely polished with AI2O3 powder, yielding a single-peaked mosaic 
structure of ~ 0.02° full width at half maximum (FWHM). The x-ray diffraction mea- 
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surements were performed on the XRD2 beamline, placed after a dipolar source at the 
Laboratorio Nacional de Luz Sincrotron, Campinas, Brazil.— The sample was mounted on 
the cold finger of a commercial closed-cycle He cryostat with a cylindrical Be window. The 
cryostat was fixed onto the Eulerian cradle of a commercial 4+2 circle diffractometer, ap- 
propriate for single crystal x-ray diffraction studies. The energy of the incident photons was 
selected by a double-bounce Si(lll) monochromator, with water-refrigeration in the first 
crystal, while the second crystal was bent for sagittal focusing. The beam was vertically 
focused or collimated by a bent Rh-coated mirror placed before the monochromator, which 
also provided filtering of high-energy photons (third and higher order harmonics). Unless 
otherwise noted, a vertically focused beam was used in our measurements, delivering, at 
7.24 keV, a flux of 3 • 10 10 photons/s at 100 mA in a spot of ~ 0.6 mm (vertical) x 2.0 
mm (horizontal) at the sample, with an energy resolution of ~ 5 eV. Our experiments were 
performed in the vertical scattering plane, i.e., perpendicular to the linear polarization of 
the incident photons. In most measurements, a solid state detector was used, except in the 
polarization study, where a scintillation detector was placed after a Ge(lll/333) analyzer 
crystal. At the energy corresponding to the Gd Ljj edge, the analyzer placed at the Ge(333) 
reflection selects (o — > a') scattering from the sample (i.e., scattered photons with the same 
polarization as the incident photons), while Ge(lll) does not significantly discriminate the 
photon polarization ((er — > a') + (a — > n') channel). 

III. MAGNETICALLY ORDERED PHASE 
A. Polarization and Resonance Properties 

Above ~ 39 K, all the observed Bragg peaks were consistent with tetragonal symmetry 
of GdRhlns (space group P4/mmm), without any detectable magnetic contribution. Below 
T/v = 39 K, additional (h, k, I) Bragg reflections (h integer; k, I half-integers) were observed. 
Such reflections were dramatically enhanced at the Gd L n and L in edges (E = 7.93 and 
7.24 keV, respectively) due to resonance phenomena (see below). Polarization analysis at the 
Gd Ln edge using a Ge(lll/333) analyser crystal demonstrated pure (a — > it') scattering at 
these fractional positions at the reciprocal space, showing that such reflections are magnetic 
in origin, with dipolar resonances at this edge . 2 ^ 2 ? 
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FIG. 1: Energy-dependence of the integrated intensity of the (0,2, |) Bragg reflection (symbols) 
across the Gd Lm and Ljj edges. Data were corrected by the absorption coefficient /i obtained 
from the fluorescence yield (solid lines). 



The energy-dependence of the absorption-corrected intensity of the (0,|, |) magnetic 
Bragg reflections at ~ 12 K are shown in Fig. Q around the Gd Ljjj and Lu edges (filled 
symbols). The energy-dependencies of the absorption coefficient, fi(E), obtained from fluo- 
rescence emission, are also given as solid lines. Resonant enhancements of over three orders 
of magnitude were observed at both edges. The intensity maximums occur ~ 2 eV above 
the absorption edges, which were defined as the inflection points of fi(E). This result is 
consistent with a dominant dipolar nature (2p — > 5d) for both resonances. Intensity oscilla- 
tions of the (0,|, |) magnetic peak were also observed above the edges, which we ascribe to 
a magnetic diffraction anomalous fine structure (DAFS). 28 
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TABLE I: Comparison between observed and calculated intensities of magnetic Bragg reflections 
at 12 K, normalized by the most intense reflection, assuming the moments fn along each one of the 
three axis of the unit cell. Experimental data were taken on resonance conditions, 2 eV above the 
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B. Magnetic structure 

The positions in reciprocal space where magnetic Bragg reflections were observed lead 
to AFM structures with propagation vector f — (0, |, |). Since there is a single magnetic 
ion per chemical unit cell in the tetragonal structure of GdRhlns, the relative neighboring 
spin orientations are unequivocally determined from the r-vector. According to this, linear 
ferromagnetic chains along the a-direction are antiferromagnetically-coupled along the other 
two axes (i.e., a C-type AFM structure). The direction of the magnetic moments may be 
obtained from the intensities of some magnetic Bragg peaks. The expression for such inten- 
sities shows a simple form in the present case for dipolar resonances and colinear magnetic 
structures, and is given by I M (f) cx (m ■ k s ) 2 , where m is the magnetic moment, and k s 
is the wave vector of the scattered light , 2 ^ 2 7 Comparisons of the observed intensities with 
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FIG. 2: Magnetic structure of GdRhlns. 



calculated ones for the magnetic moments along a, b, and c-directions are given in Table Ufor 
both compounds. Good agreement between calculated and experimental data are obtained 
for the magnetic moments pointing towards the direction of the ferromagnetic chains, i.e., 
the a-direction. Figure 121 displays the magnetic structure of GdRhIn 5 , as determined above. 
This structure is consistent with a 155 Gd Mosbauer spectroscopy study in this compound, 
which indicated a colinear magnetic structure with the spins lying in the ab plane.— 

C. Search for symmetry lowering of the crystal or electronic structure 

It is interesting to note that the magnetic structure shown in Fig. |2] does not have 
a tetragonal symmetry, in the sense that the ferromagnetic chains are aligned along one 
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TABLE II: Extraction of the a, b, and c lattice parameters at 20 K and 300 K from selected 
(00Z), (hOl), and (Okl) Bragg peaks. The width of radial (9 — 29) scans are also given in terms of 
59 /tan{6) (full width at half maximum), were 9 is the Bragg angle. 
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specific direction in the ab plane, namely the a-direction. It is therefore valid to ask whether 
a symmetry lowering of the crystal and/or electronic structure from the P4/mmm space 
group also occurs, either by an orthorhombic distortion of the lattice parameters or by the 
presence of a charge density wave. These possibilities were also investigated by synchrotron 
x-ray diffraction. 

Firstly, a hypothetical orthorhombic distortion was probed by measuring the Bragg inter- 
plane distances dhki of a set of (00/), (hOl), and (Okl) charge reflections, using a monochro- 
matic x-ray beam with energy E = 14472 eV and resolution (5E/E = 1.0 x 10~ 3 ). The 
irradiated (001) surface was the same used in the magnetic diffraction investigation. A scin- 
tillation detector was placed after a Si(lll) analyzer crystal. We used the usual relation 
d\ kl = l/[(h/a) 2 + (k/b) 2 + (l/c) 2 ] for orthogonal axes, obtaining the lattice parameter c 
directly from <i o«- The a and b parameters were then obtained from dhoi and d^i- Table 
HTl shows the results at 20 K and 300 K, indicating a metrically tetragonal phase within 
experimental errors. 

The possibility of an orthorhombic structure presenting a mosaic of domains with inter- 
changed a and b axes was also considered. In this case, a two-peak structure is expected 
at each (hOl) or (Okl) reflection. Such a structure was not observed in our measurements, 
rather yielding single symmetric Lorentzian line shapes. In this case, the upper limit for the 
orthorhombic distortion is set by the measured peak widths. This information is also given 
in Table ILTl in terms of S9/tan(9). A larger S9/tan(9) for (hOl) and (Okl) reflections with re- 
spect to (00/) reflections would be consistent with an orthorhombic distortion or anisotropic 
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strain in the ab plane. Nonetheless, it is seen that no additional broadening of the (hOl) or 
(Okl) peaks with respect to (00/) was observed within our resolution. In addition, the peak 
widths are the same at 20 K and 300 K, showing the absence of any temperature-dependent 
distortion or strain. Using the data of Table El the upper limit of the hypothetical or- 
thorhombic distortion is inferred to be ((b — a) /a) < 2 x 10~ 4 for both single-domain or 
mosaic (multi-domain) distortions. We conclude that no relevant bulk lattice distortion or 
anisotropic strain associated with the anisotropy of the magnetic structure along the a-axis 
take place in this compound. 

The possible presence of charge density waves (CDWs) was probed by a systematic search 
in reciprocal space at x-ray energies of 7930 eV (resonant condition at Gd Lu absorption 
edge) and at 7106.7 eV (below the Gd L IH absorption edge). A set of charge reflections 
- (001), (002), (011), (022), (021), (102) - was chosen to define the borders of the one- 
dimensional scans in reciprocal space. Several scans along high symmetry h, k, and I mixed 
and unmixed directions were performed at T = 11.5 K. The search was completed by two- 
dimensional h — k, h — I and k — I maps with a common border at the (022) reciprocal 
space position (Sh,k,l = 0.5). We did not find any extra peak that might be assigned to 
previously unknown ordered structures. Although our scan procedure had been unable to 
directly reveal the existence of new electronic or structural phases in GdRhlns, we cannot 
completely rule out the possibility of existence of low-symmetry incommensurate phases 
such as CDWs, since we did not investigate the whole reciprocal space volume. 

D. Discussion: Rise of the C-AFM magnetic structure and competition with G- 
AFM: the role of long-range exchange interactions 

In the i? m (Co,Rh,Ir) n In 3m+2n (R = rare earth ^ Ce) family, the main coupling between 
the R ions comes from the Ruderman-Kittel-Kasuya-Yosida (RKKY) mechanism.— If we 
consider antiferromagnetic couplings only between the first- (Ji) and the second- nearest- 
neighbors (J2), the magnetic ground state is determined by their relative strength a = J\j J 2 . 
The G-AFM structure represented in Figure El frustrates J 2 but satisfies J\ and would be 
favored by a >> 1. On the other hand, the C-AFM structure satisfies part of J\ and 
all J 2 interactions and thus may be favored if J\ and J 2 are comparable in magnitude. 
For the RKKY interaction, a usually depends on the first- and second-neighbor distances 
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C-AFM G-AFM 

FIG. 3: (Color online) Two-dimensional representation of the C-AFM and G-AFM spin structures. 
The first- and second-nearest-neighbor exchange interactions (Ji and J 2 , respectively) are also 
indicated. 



among the rare-earth ions and on the topology of the Fermi surface of each compound. 
However, in these compounds the topology and volume of the Fermi surface depends little 
on the particular rare-earth ion since, except for the Ce based compounds, the localized 4/ 
electrons do not hybridize considerably with the conduction band. Also, we note that the 
tetragonal a-lattice parameter shows only small variations along the -R m (Co,Rh,Ir) n In3 m+2ri 
series (< 2 %).^£ It is therefore not completely surprising that the magnetic structures of 
all known i? m (Co,Rh,Ir) n In 3m+2n (R = Nd to Gd) compounds are similar, characterized 
by ferromagnetic chains along a specific first-nearest-neighbor direction (a) with an AFM 
coupling along b and c, i.e., the C-AFM structure (see Fig. Eland refs . 19 ' 20 ' 21 ' 22 ). 

The orientation of the magnetic moment in any of these magnetic structures is deter- 
mined by anisotropic interactions. The direction of the staggered moment on GdRhIn 5 is in 
agreement with other tetragonal compounds with the same antiferromagnetic wave-vector 
and staggered moment direction, such as GdAu^i^ and GdCu 2 Si 2 .— It has been argued 
that in theses cases the magnetic dipolar interaction is the dominant source of anisotropy 
(being of the order of tens of fieV)M The same appear to hold for GdRhlns. Nonetheless, we 
should mention that other possible sources of anisotropy for Gd compounds, such as a spin- 
orbit coupling of the conduction electrons 3 ^ or crystal electric field via excited state d 35 ' 36 ' 37 ' 38 
might in principle be relevant to the problem. 
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The same general conclusions for the ground state may remain valid when the spin 
anisotropy due to the relevant crystal field effects is introduced in the Hamiltonian for 
R 7^ Gd. In this case, the crystal field effects determine the spin direction with respect 
to the unit cell axes and may affect T^,— but the relative orientation between neighboring 
i?-spins is still determined mainly by a. 

The extension of the above scenario for the Ce-based compounds is not straightforward. 
This is because the Ce 4/ electrons may be hybridized with the conduction band. In cases 
where the Ce 4/ electrons are itinerant, the RKKY mechanism is no longer applicable. On 
the other hand, for compounds with localized 4/ moments, the above scenario of J\j Ji com- 
petition might be useful. Particularly, deHaas van Alphen measurements on Cei_ x La x RhIn 5 
showed no significant change in the Fermi surface topology or volume over the entire doping 
range (0 < x < l), 40 showing that the 4/ electrons remain localized. Thus, we focus our 
discussion on the magnetic structure of CeRhIn 5 with localized moments. This compound 
shows AFM alignment along two directions (a and b), and a spiral alignment along the 
tetragonal c-axis,— defining a propagation vector f = (\, \, 0.297). This magnetic structure 
may be seen as an intermediate case between the G-AFM structure (f[ = (|, |, |)) and a 
C-AFM structure with the FM chains along the 5-direction (f 2 = (|, \ , 0)). Notice that the 
C-AFM structure is expected to be a competitive ground state, since a is expected to be 
similar for CeRhIn 5 and GdRhIn 5 due to presumably similar Fermi surfaces. On the other 
hand, the G-AFM state also appears to be competitive for the Ce-based compounds, since it 
is the ground state of Cdn a 1 ^ 11 and C^Rhlng.™ It is therefore not implausive to infer that 
the incommensurate magnetic structure of CeRhlns is actually a result of a close competi- 
tion between the G-AFM and C-AFM ground states, paving the way for the stabilization 
of an intermediate spiral phase, perhaps with the aid of very long-range RKKY interactions 
(J 3 , J 4 , etc). 

IV. CRITICAL BEHAVIOR 
A. Experimental 

The critical behavior of the sublattice magnetization at the paramagnetic transition was 
investigated. FiguresEfa,b) show the T-dependence of the intensity of the (0, |, y) magnetic 
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FIG. 4: (Color online) Temperature-dependence of the (a,b) integrated intensity and (c) width of 
radial (6 — 29) scans of the (0, \, -4-) magnetic Bragg reflection of GdRhlns and (1,0,4) reflection 
of Gd2lrlng (taken from refi^2). A fit to a critical power law behavior, (T/v — T) 2 ^, characteristic 
of a second-order transition, is given in (b) as a line. The inset show the correlation length L, 
obtained from the data in (c) after deconvolution of the instrumental width. The experiment was 
performed in resonance conditions, 2 eV above the Gd L//-edge. 

Bragg peak of GdRhIn 5 in different T-intervals. Close to and below Tv = 39 K, this 
could be fitted by a power- law behavior, / oc (1 — T/T^) 213 , which is characteristic of a 
second-order transition. The experimental determination of the critical parameter (3 depends 
slightly on the temperature interval in which the fitting is performed. For fits in T-intervals 
(Tv — T) /Tv < 0.01, < 0.03, and < 0.05, one obtains (3 = 0.370(16), 0.346(5), and 0.339(4), 
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respectively. Quoted errors in parentheses are statistical only, and represent one standard 
deviation. Figure Efc) shows the T-dependence of the width of the (0, |, ^) magnetic peak, 
as obtained in radial (9 — 28) scans. Much below T/v, the width is instrumental only, 
indicating long-range order with correlation length above ~ 1000 A. For 0.998T/V < T < 
T/v, a peak broadening above the instrumental resolution was noticed. The inset of Fig. 
Efc) displays the magnetic correlation length in this T-interval, obtained from the peak 
broadening data. We should mention that, above T/v, the intensities were below our detection 
limit for this sample, thus the short-range dynamic correlations in the paramagnetic phase 
could not be investigated. 

These results for GdRhlns may be compared to the previously reported measurements for 
Gd2irlng, where a larger critical parameter f3 = 0.39 was obtained for the critical interval 
(T/v — T)/T/v < 0.10.— The comparison becomes clear in Fig. EJb), where the critical 
behavior of Gd2irln 8 , using data of ref.— , is directly compared to GdRhIn 5 . In the ordered 
phase (T < T/v), both compounds appear to show identical behavior, described by the 
same critical exponent (3 ~ 0.35. However, while the magnetic intensities tends to zero as 
T ~ T/v for GdRhlns, a significant residual scattering was observed near and above the 
transition for Gd2irln 8 , which may be mostly ascribed to magnetic correlations in the near 
surface region due to long-range correlated quenched disorder . 22 i 41 This effect smooths out 
the transition observed by x-ray diffraction, and interfere severely in the extraction of the f3 
critical exponent. We conclude that, although both compounds were equally finely polished 
before measurements, our studied Gd2irln 8 surface showed a larger degree of near-surface 
disorder, leading to a less reliable extraction of the f3 exponent. Perhaps the most relevant 
information from this analysis is that a direct comparison of the magnetic intensities of 
Gd2irin 8 and Gd2RhIn 5 near T/v reveals an identical behavior in the T-interval where the 
surface disorder effects are negligible, consistent with an identical critical exponent j3 ~ 0.35. 
This conclusion is consistent with the prediction that both compounds belong to the same 
universality class for magnetism (see below). 

B. Renormalization Group Analysis 

The experimental critical behavior described above may be compared with theoretical 
expectations based on the symmetry of the crystal and magnetic structures of this compound. 
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As demonstrated by Mukamel et a/. , 42 > 43 the Ginzburg-Landau- Wilson (GLW) Hamiltonian 
for the class of antiferromagnetic problems is written in terms of a staggered order parameter 
(OP) with a total number of components N = nm, where n is the number of spin components 
allowed by the irreducible representations of the paramagnetic space group, or equivalently, 
the number of degrees of freedom for the critical spin fluctuations, and m is the number of 
ways the unitary cell may be enlarged by all distinct AFM ordering wave-vectors allowed by 
the symmetry of the crystal. The total number of components may easily exceed iV > 4, 44 
opening up the possibility to classical fluctuation-induced first-order transitions, such as 
those already reported experimentally.—" 4 ^ The criterion for the study of phase transitions 
using renormalization group (RG) theoretical methods is based on the stability of the fixed 
points in the RG flow, in the sense that the phase transition is of second order when the 
flow continuously approaches a fixed point which is stable with respect to the fluctuations 
of the staggered field, and indicates something else like an abrupt or an smeared transition 
otherwise, when the flow exhibits a runaway.— 

Because of the critical exponents universality, the phase transition is sensitive only to a 
few parameters like the dimension d of the system, the number of components N of the OP 
and the symmetry of the crystal, which reflects in the anisotropies of the Hamiltonian. As 
shown by Brezin et aZ.jsa for iV > 4, the isotropic fixed point is always unstable with respect 
to anisotropies in the original Hamiltonian, while for N < A this fixed point is always stable 
because the O(N) symmetry is dynamically generated near the critical point. 

Here, we concentrate the analysis on a few compounds of the Gd m (Rh,Ir) r Jn 3m+ 2n series, 
including the title compound. In the cubic Gdln3 crystal, the observed AFM propagation 
vector t = |, 0)^ represents an ordering state which is degenerated with the states 
represented by the propagation vectors (0, |, |) and (|,0, |), resulting in m = 3 distinct 
AFM propagation vectors allowed by the crystal symmetry. Since n < 3 is the total number 
of spin degrees of freedom, the staggered OP of Gdln3 has a total number of N = 3n < 9 
components. In the tetragonal Gd2lrln§ crystal, the cubic symmetry is broken in the [001] 
direction by non-magnetic planes of Ir. 22 This crystal orders antiferromagnetically along the 
t = (|,0, 0) direction, which is equivalent by symmetry to the (0, |,0) direction, giving 
m = 2 and N = 2n < 6 OP components. Since the magnetization of GdRhIn 5 orders along 
the t = (0, |, |) direction (which is equivalent by symmetry to the (|,0, |) direction), its 
staggered OP has the same number of components of GdIr 2 In 8 , N = 2n < 6. 
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In isotropic magnetic systems, the equality n = 3 holds for the number of spin degrees 
of freedom, and all above relations for N are valid with the equality sign. We conclude 
that if GdRhIn 5 , Gd2lrln 8 , and Gdln 3 were perfectly isotropic spin systems, first-order 
transitions would be obtained for all three compounds,— in contrast with the second-order 
transitions observed for GdRhIn 5 (see Fig. HJ) and Gd 2 IrIn 8 (see ref.— ). We conclude that 
spin anisotropy must be properly taken into account for a correct analysis of the critical 
behavior or these compounds. As noted in Section III.D, a number of distinct possible 
sources of anisotropy may be anticipated for Gd compounds with half-filled 4/ 7 shell, such 
as dipolar interactions,— spin-orbit coupling of the conduction electrons^ or crystal electric 
field via excited states . 35 i 36 i 37 i 38 The strength of such interactions is typically of the order of 
tens of /xeV . 32 i 34 Even though this energy scale is about three orders of magnitude weaker 
than typical exchange energies in Gd systems, we note, for example, that there are consistent 
evidences that in Gd metal the dipolar interaction is not only responsible for the ground 
state anisotropy but also determines its critical behavior . > 50 

In order to proceed with our analysis, the major source of anisotropy must be identi- 
fied. Since there are recent indications that dipolar interactions are the major source of 
anisotropy 3 ^ and responsible for the specific heat behavior in an extensive variety of Gd 
compounds , 47 i 51 i 52 we pay attention to the possible influence of the dipolar coupling in the 
critical behavior of these materials. The dipolar anisotropy breaks the rotational symmetry 
of the Gd spins by lowering the size of the space of degenerated states where the spins are 
allowed to fluctuate. In the RG sense, the influence of the dipolar interaction will be decisive 
if it proves to be a relevant source of anisotropy in a previously isotropic Hamiltonian, 53 The 
difficulty of the RG method here is that it leads to rather inconclusive results when the RG 
flow has no stable fixed points, since the flow rapidly moves towards a region where the tech- 
nique is no longer valid. To see this, we write down the most general classical Hamiltonian 
that describes the physics of the isotropic Gd spin problem, which corresponds to a GLW 
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Hamiltonian of m coupled O(n) symmetric models, 



a,i,j a^fi ij 



+W^2^2<Pai <Pf3i 4>aj <f>/3j>, ( 1 ) 

ij ) 

where a, /3 — 1, ...,m indexes the distinct AF wave-vectors and i, j = 1, ...n < 3 labels the 

spin components in a given orthogonal basis, like x, y, z. Note that all terms are written as 

powers of scalar products of the spin components <^ ■ <\) v because of the assumed rotational 

symmetry of the spins near the phase transition. The quartic terms differ only by the 

different ways to combine the Greek indexes that label the equivalent AF wave-vectors. 

Including a general dipolar interaction term, which in the antiferromagnetic case has the 
form,j£ 



Hn{<t>) = £ 



(V0 iQ ) 2 - f(di(f)i a ) 2 + h <),o„J) r ., 



(2) 



where / and h are proportional to the dipolar coupling constant (^/i^) 2 , with g[i B being the 
total magnetic moment of the Gd ion, we show in Appendix A that the total Hamiltonian 

H = H + H D (3) 

has the isotropic fixed point for iV = nm < 4, and no stable fixed points for iV > 4. In 
particular, the scaling of the dipolar parameters from equation 

% ~ 

§ - (4) 

with 

32 / l \ 2 

Vf = r]h = YVs^J ( 2 ^ + (^-i)K + ^ 2 )) >o, 

indicates that the dipolar interaction is irrelevant in the vicinity of the isotropic fixed point 
for N < 4, what corresponds in our case (n = 3) to m — 1. For larger m, however (m = 2 
for GdRhIn 5 , Gd2lrln 8 and m = 3 for Gdln 3 ), this Hamiltonian has no stable fixed points 
and, although this would point to the direction of a fluctuation induced first order phase 

17 



transition, it is not clear what happens in this case. The fact that the / and h parameters 
flow initially to zero is expected, since we do not include sources of anisotropy in the bare 
Hamiltonian as the cubic term (v) from reference^. Nevertheless, as pointed out in this 
reference, such terms could well be self generated in a higher loop expansion. 

From mean field calculations of the propagation vector and moment direction in several 
Gd compounds, Rotter et alM have shown strong evidence that in these compounds the 
observed anisotropy stems from the magnetic dipolar interaction. If we assume that dipolar 
interaction is relevant at the phase transition, it must break the spin rotational invariance 
of Hamiltonian (JI}. Once this symmetry is broken, the quartic terms should be written 
in the most general way allowed by the crystallographic group. In the particular case of 
the Gd series compounds we have studied, we have not identified bi-critical, tri-critical or 
multi-critical points associated to the paramagnetic phase transition, meaning that the spin 
fluctuations are confined inside subspaces of degenerated spin configurations. In other words, 
each of these subspaces in the spin space correspond to an irreducible representation of the 
AF order parameter. The number and the size of all the irreducible representations allowed 
follows directly from the crystallographic point group symmetry and from the position of the 
ordering wave-vector f in the Brillouin zone (BZ). In the case of GdRhIn 5 , the P4/mmm 
space group associated with the special point r = (|,0, |), at the border of the BZ, has 
three irreducible representations of dimension m, with the spin pointing along the unit-cell 
axes of the crystal (for details, see Appendix B). In this case, the number of spin degrees of 
freedom for each representation is just n = 1. The most general Hamiltonian would be 

H = \ £ M« + c W«) 2 ] + u £ <t>l€ 

a a 

+«£#$. (5) 

with the dipolar interaction in this case simply renormalizing the gradient term. For m — 2, 
as in the case of GdRhIn 5 , this Hamiltonian has one 0(m) symmetric fixed point^ with 
(3 = 0.36 in two-loop expansion. The ordering wave vector r = (i,0,0) of Gd2irln 8 (also 
from space group P4/mmm) gives the same three irreducible representations of dimension 
m, and therefore GdRhIn 5 and Gd2irln 8 both lie in the same universality class. In Gdln 3 , 
the Pm3m group acting on a BZ with the ordering wave-vector f = (~, |,0) admits one 
irreducible representation of size N = m, with the spin pointing along the c-axis, and 
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another one of size N = 2m with the spins confined in the ab plane. The Hamiltonian of 
first representation follows the general form of Eq. (J3J, showing one isotropic fixed point 
for m < 4. The second representation is larger (N = 6) and its Hamiltonian follows Eq. 
(J3J), where no stable fixed points were found, which would indicate the possibility of a 
fluctuation-induced first-order transition for this case, even if the anisotropy is included. 

V. CONCLUSIONS 

In summary, our resonant x-ray diffraction experiments show a C-AFM magnetic struc- 
ture for GdRhlns with the spins lying along the FM chain direction. This structure is 
rationalized in terms of a competition between first- and second-neighbors exchange inter- 
actions. This scenario was extended to other members of the i? m (Co,Rh,Ir) n In3 m+ 2 ra (R = 
rare earth) family, in particular the Ce-based heavy-fermion superconductors. The critical 
behavior close to the paramagnetic transition was investigated, revealing a second-order 
transition with critical exponent (3 ~ 0.35 for GdRhlns and Gd2lrlng. A renormalization 
group analysis predicts that both compounds belong to the same universality class, in agree- 
ment with experiment. However, a first-order transition is predicted in the absence of spin 
anisotropy terms in the Hamiltonian, in contrast to our results, indicating that such inter- 
actions may be important to stabilize a critical point in this family. Indeed, we show that 
if dipolar interactions or any other relevant source of spin anisotropy that allows the Gd 
spins to couple with the lattice is included, then the renormalization group analysis pre- 
dicts a second order phase transition for GdRhIn 5 and GdIr 2 In 8 with f3 = 0.36 in two-loop 
e-expansion, in good agreement with our experiments. 

APPENDIX A: RG PROCEDURE 

Since the dipolar interaction is marginally relevant at the one loop level, the RG calcula- 
tion has to go to a second loop expansion. Following the standard RG procedure along the 
lines of reference,^ we integrate out of the partition function the fluctuation modes with 
wave-vectors b^ 1 < |k| < 1 (which we denote by J dk and 6 > 1) and rescale the spins. 
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This leads to a renormalization of the dipolar part of the Hamiltonian (J2J) according to 

T 2 (q) = -16 (2m 2 + (m - 1) (v 2 + w 2 )) x 
y dkGjj (k) ^2 hmm (k - q) 

-32 (2u 2 + (m - x 

/ dk£G 4m (k)£j lmi (k-q) (Al) 



where 



/ 7 */j (p) = / d d kG 75 (k) G 7/3 (k + p) , 

y&— i <iici<i 



'6- 1 <|k|<l 

We assume that the dipolar interaction is strong enough in the critical region, i.e T — T c <C 
^ gM f^ , so that the bare propagator is given by 



G l3 (k) = 1 



^■-Vjj/ + /o ( ^ ) ^ 



with /o and /i-o being proportional to the dipolar coupling constant and to geometric factors 
reflecting the lattice symmetry (see eq. 23 from ref.— ). After performing the integrals on 
equation (jAljl and expanding T 2 (q) to second order in q, which though rather cumbersome 
is straightforward, we obtain equation (j3J. Since / and h flow initially to zero, the available 
fixed points up to this order correspond to the fixed points of the isotropic Hamiltonian 
(0). The RG of the isotropic problem has been studied in detail by Ref.—, which indicates 
the absence of stable fixed points for = 3m > 4. The absence of stable fixed points 
within the isotropic model remains inaltered with the inclusion of long-range correlated 
quenched disorder.^ In one-loop at least, the disordered isotropic model produces no new 
stable fixed points 5 - beyond the two unphysical fixed points previously found by Halpering 
and Weinrib. 1 *^ 



APPENDIX B: IRREDUCIBLE REPRESENTATIONS OF GdRhIn 5 

The simplest way to obtain all the irreducible representations allowed by the crystallo- 
graphic point group for a given ordering point r in the BZ is to decompose the spins in a 
given basis fixed with respect to f, and then apply all the point group symmetry operations 
to see how the spin components in the original basis will change. The idea is that if we 
properly choose the original spin basis, then we are able to identify the subspaces where the 
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spins will be confined by the application of the symmetry operations allowed by the crystal 
point group only. For GdRhlns, which ordering wave- vectors are 

.1 1. 

T«=l = (3,0,-) 
T a=2 = (0,-,-) 

we will define the spin basis by #2,a#3,a}, where 

u l>a=1 = (1,0,0) 

#2,0=1 = (0,1,0) 

#3,0=1 = (0,0,1) 



Z?l,a=2 = (0,1,0) 
#2,a=2 = (-1,0,0) 
#3,a=2 = (0,0,1). 

Denoting 0j jCe as 0^1 = 0, and 0j.2 = 0i for the i-th spin component with respect to the a-th 
ordering wave vector, the symmetry operations of the P4/mmm point group generators are 

C 4 [001] : fa 
C 2 [100] : 0! 

01 

i : (Pj 

We see that once a spin points along one of the principal axes (labeled by the i index) the 
application of the crystallographic point group symmetry operations will "trap" it on the 
same direction. This results in 3 irreducible representations of size N = m = 2. 



«-> 01, 02 <-> 02, 03 ► 03 ► -03, 

> -01, 02 > 02, 03 > 03, 

► -01, 02 > 02, 03 ► -03, 

— > -0i> 0j — -0j- 
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